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; Abstract 
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We show that the phenomenological AC 2 M 2 ansatz is consistent with QCD through 
order 1/rrib in the description of B —>■ Ivi + X u and B — > 7 + X s transitions, including 
their energy spectra and differential distributions. This suggests a concrete realiza- 
tion for the QCD distribution function, which we call the "Roman" function. On the 
other hand the AC 2 M 2 model description of the end-point domain in B —>■ lv\ + X c 
is incompatible with QCD: a different distribution function enters the description of 
b — > c decays as compared to the transitions to the massless quarks. Both observa- 
tions - the validity of the AC 2 M 2 -like description for heavy-to-light transitions and 
the emergence of the new distribution function in the b —>■ c case - are in contradiction 
to a recent claim in the literature. The intrinsic limitation of the AC 2 M 2 model could 
reveal itself in different values of the effective b quark mass from fits of the A^ and B 
decays. 
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1. In recent publications [[T], H we have shown how the lepton and the photon 
energy spectra in semileptonic and radiative decays, respectively, of beauty hadrons 
can be predicted from QCD proper without recourse to phenomenological ad hoc 
assumptions; this is achieved by employing a heavy quark expansion in inverse powers 
of the beauty quark mass to deal with non-perturbative corrections. It turned out 
that the shape of the spectrum in the end-point region requires a special treatment 
(see also Refs. [Q, [|). The physical phenomenon behind that is quite transparent: 
the heavy quark Q moves inside the decaying hadron Hq and this motion smears 
the lepton spectrum, most significantly in the end-point region. It has been noted 
|], that a rigorous QCD description - though having a similar origin - does not 
generically coincide with simple non-relativistic models of the 'Fermi' motion where 
the potential energy of the light degrees of freedom is always small as compared to 
their space-like momentum; the most popular model of this type is due to AC 2 M 2 ||. 
Original generalities of the QCD approach are given in Refs. 0, 0, § • 

The AC 2 M 2 model describes the beauty meson B in very simple terms. It treats it 
as consisting of the heavy b quark plus a spectator with fixed mass m sp ; the latter is 
of the order of a typical hadronic scale and represents a fit parameter. The spectator 
quark has a momentum distribution where p is its three-dimensional momen- 

tum. The b quark then cannot possess a fixed mass. Indeed, the energy-momentum 
conservation implies for such a simple boundstate that the b quark energy Ej, is given 
by Mb — \J\p\ 2 + ^ 2 P - The mass of the b quark equals E b through order p; thus one 
obtains a "floating" b mass 




which depends on \p\. The lepton spectrum is first obtained from b quark decay in 
a moving frame where its momentum is — p\ the quark mass is assumed to be m(. 
Then the distribution is averaged with the weight function (further details can 

be found in the original work ||). The Galilean transformation above introduces a 
Doppler smearing of the spectrum that is of order l/m^ the transverse Doppler shift 
arises only at order 1/m 2 and is ignored here. 

The AC 2 M 2 model is extensively used in the analysis of the lepton energy spectrum 
in semi-leptonic decays. Surprising though it is, given all the naivete of the model, it 
reproduces very well numerically the shape of the semi-leptonic spectra in its regular 
part, as it is derived from QCD (see In this paper we analyze in more detail 
what features of QCD are reflected in this model for the end-point spectrum. We 
will show that - contrary to a first impression - this model is built sufficiently well 
to reproduce correctly the pattern of the QCD description of the heavy quark motion 
provided one considers the transition of the heavy to a massless quark. Moreover, 
the specific choice adopted in AC 2 M 2 for the weight function $(|p|) corresponds 
to a particular form of the genuine distribution function F(x) introduced in Ref. 
which schematically can be viewed as 

4 *» 1 
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where p + is the light cone combination po + p z entering the QCD description (the 
exact relation to F(x) will be stated below). This represents an acceptable realization 
of the QCD distribution function meeting all necessary requirements (we will call it 
the Roman function []) ; it can easily accommodate improvements should they become 
necessary. At the same time the AC 2 M 2 model per se does not reproduce correctly 
the dependence on the final state quark mass m q if the mass is large enough. As 
follows from the QCD analysis M the boundary lies at m q ~ (Amj,) 1 / 2 and c quark is 
sufficiently heavy in this sense. Moreover, in the b — > c transitions the AC 2 M 2 model 
in its original form would yield spectra beyond the physical kinematic boundary. A 
good fit of the experimental lepton energy spectrum in b —>■ c is achieved [|J at a price 
of introducing an upper cut off in the heavy quark momentum depending on the final 
quark mass. 

Recently a comparison of the QCD result and of the AC 2 M 2 model has been given 
in Refs. || |i~0fl , where it was pointed out that the lepton energy spectrum obtained 
from the AC 2 M 2 ansatz can be considered to be consistent with the QCD prediction 
in the sense that there are no corrections to order l/m^. In this note we give a more 
detailed analysis. We also point out that the consideration of Refs. || f|] in the part 
concerning the AC 2 M 2 model is incorrect. 

2. Some of the dynamical features in the spectra of charged leptons get obscured 
by the integration over the neutrino energy. They are seen in a cleaner way in the 
spectrum of photons in the radiative B — > 7 + X s transitions. This spectrum in the 
free quark approximation is just a monochromatic line: 

— — — = — mld(E 3 

dE 4tt b V 2 ; v ; 

where E is the photon energy, is the (current) b quark mass, the strange quark is 
considered as massless and A is the coupling constant determining the bsj vertex as 
it emerges from the electroweak theory: 

£\b^ s ~f = -XF^sil + 75)0-^6. (4) 

In both the QCD approach [[I], 0, |J and in the AC 2 M 2 model H the monochromatic 
line of eq. (|3|) is transformed into a peak of a finite width due to the heavy quark 
motion. The width of the peak is of the order of a typical hadronic scale. Note, that 
in the AC 2 M 2 model the variation of the b quark mass (see below) also contributes 
to the width. 

Gluon radiation will also lead to a smearing out of the photon line. Since our 
primary object here is to discuss non-perturbative effects like the Fermi motion, we 
will ignore perturbative gluon emission throughout most of this paper. 

Non-perturbative corrections are incorporated in the AC 2 M 2 model through the 
three fit parameters m sp , m q and pp. Let us recall now some general expressions 

1 The possible shape of the function shown in Figs. 1,2 demonstrates that it is definitely not of a sans serif 
type. 
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derived previously from QCD for m q = case. Neglecting the gluon radiative 
corrections one can write the spectrum for the inclusive radiative decay in the form 

where A = Mb — nib, 

x = l(2E-m b ) , (6) 
and F(x) is the QCD distribution function. The moments of F(x) , 

dxx n F(x) , n = 0,l,... (7) 

are related to the expectation values of local operators taken between B mesons: 

(B\ S Im^-'K^b- traces \B) = a n k n {v^...v^ n - traces). (8) 



2M F 



Here S is the symmetrization symbol and t> M is the four-velocity of the B meson. 
Moreover, 

tt m = iDf, - rribV^ (9) 

where _D M is the covariant derivative. Note that ao = 1 and a\ — (these relations hold 
up to 1/m 2 corrections which we disregard throughout). An important consequence 
of quantum chromodynamics is that the total width 

r = ^mg + ... (10) 

receives no corrections at the level l/m b ; they arise first at order 1/m 2 |7|, §. 

It should be noted that the true distribution function F(x) is always one-dimensio- 
nal. For in the QCD description different components of the momentum represented 
by 7r^ do not commute and the distribution function can be introduced only for a 
single combination of the components of the momentum. In the case at hand, when 
the final quark is massless, i.e. m q = 0, it is the 'light cone' distribution function 
that enters (i.e. we actually deal with the distribution in ttq + tt z where the z-axis is 
chosen along the direction of the photon momentum; for more details see Ref. 0). 

At first sight the two approaches - that of Refs. [I], 0, |J on the one hand, and 
the AC 2 M 2 model on the other - have completely different properties. Two main 
distinctions are obvious. 

The AC 2 M 2 model introduces the distribution over the three-dimensional nonrel- 
ativistic momentum, $(p) where p = \p\ is the absolute value of the b momentum. 
This distribution results in the Doppler smearing of the photon spectrum with the 
spread linear in l/mt,. Moreover, the AC 2 M 2 prescription replaces m;, by m(, the 
floating quark mass, see eq. (|l|): 



dT ^M A 2 f a , _ ^ a 3 . m P z^ 

~JE~ = 8^J P P ^ ( m b) 6 ( E ~2~Y^' ( ^ 
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where z = cos 6 and 9 is the angle between p and the direction of the photon momen- 
tum. (Eq. ( |XT| ) is obtained by substituting m& by m{ , making the Lorentz boost and 
then co nvo luting the spectrum resulting from the parton approximation with 

At the same time, the QCD expansion operates in terms of the fixed, current 
mass of the the b quark, m&. The effect of the Fermi motion enters via the light cone 
distribution function defined in eqs. (|7|), (^) (we remind the reader that the final s 
quark is treated as massless like the q quark before). 

3. We will show now that in spite of these differences the two approaches can 
actually lead to the same photon spectrum. The AC 2 M 2 ansatz yields for the photon 
spectrum after averaging over the direction of p 



acm A 2 



0(e)/ dp 2 §{p 2 )(M B -Jm% + pA , (12) 




dE 8n J P 2 (e) 

with 

e = Mb — 2E , p*(e) = -\^- e j . (13) 

The spectrum is automatically cut off at the true kinematic boundary, i.e. at E — 
Mb/2 by the step function for e < 0, irrespective of the particular choice of $(p). 
This is not surprising, of course, because the adoption of the floating b quark mass 
was dictated just by this requirement. 

The shape of the spectrum is obtained by direct integration and depends on the 
specific form of $(p). The distribution $(p) suggested in Ref. and routinely used 
in experimental analyses is 

4 -4 

$(p) = ^^e (14) 
Pfv 71 

with the following normalization 

/•oo 

dpp 2 &(p) = 1. 







Notice that the upper limit of integration (which must be Mb) is extended to oo. This 
introduces only an exponentially small error; below we will always use this extension 
of the upper limit of integration. 

Using the ansatz of eq. (0) and keeping terms through order pp/Ms one gets: 



dT ACM (X 2 M B S 



0(e) 



P 2 (e) 



dE \ 4tt J Pf\[t^ 
where the energy E enters via e, pl(e) defined in eq. fljjD and 



M b ' { 2' P+ p| ,+U[ Mj 



B 



(15) 



2 

m 



p=-f. (16) 
Pf 

The factor in the square brackets can be put to unity in discussing the shape of 
the spectrum in this approximation. The term proportional to pf/Mb includes the 
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incomplete gamma function T(a,x) and will be taken into account below to obtain 
the total width to l/m& accuracy. The expression in eq. (|15|) exhibits a pronounced 
peak whose shape depends, of course, on the value of p (see Fig. 1). For p < 1 it 
is rather asymmetric, a second gratifying feature of the AC 2 M 2 model, since it is in 
qualitative accord with the findings in QCD. As we will see later on, this model can 
provide an approximate description of the real world only when p is rather small. 
Straightforward integration of the spectrum eq. (JT2|) leads to 



\ 2 Ml 

ACM — 



ZPFP_ , 2 p ) + 0( 1& 



(17) 



where K\ is the McDonald function. Once Tacm is expressed in terms of the quantity 

m b 



^CM = Mb _ PFP^^ 



7T V 2' 



(which, as we will see in a moment, is nothing but the value of the floating mass m[ 
of eq. ([!]) averaged over the distribution $(p) ), the correction to first order in l/ra& 
is eliminated from the total width, 



X 2 (m 

' ' ACM = ~ 



ACM\3 



47T 



l + 0(l/m 2 )), (19) 



in full agreement with the general statement of the absence of the l/m& correction in 
the total width [Jl IB, EL Thus, one must identify 

M B - mb ^^K4^[^ at^. (20) 



That the first order correction to the total width can be absorbed into model 
parameters is not surprising. What is less obvious is that this definition is compatible 
with QCD predictions for other quantities as well. We will show now that this is 
indeed the case. (The same conclusion concerning inclusive semi-leptonic decays has 
been reached previously in Refs. || |i0|.) 



From the QCD analysis of Refs. |2], |TT[ we know the average photon energy that 
points to the true mass m b : 

(e) = (M B — 2E) = A(l — (x)) = A (21) 

where x is defined by eq. (^ and (x) = a± = 0. To find (e) in the AC 2 M 2 model we 
use eq. ( |12"D or eq. ([15]). Now the difference between the quark mass and Mb can be 
ignored and one obtains 

(e) ACM = ^e^K^) (1 + 0{l/m b )) . (22) 

'7T 2 



which agrees with the QCD expression eq. (^T|) provided that the identification of 
eq. (p0| ) is made. Moreover, this holds true for any choice of 
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For the second moment of the photon energy one derives from eq. (12) 

(e 2 )ACM=p 2 F (2 + p)- (23) 
On the other hand, according to Ref. 

<e 2 >=A 2 (l + <* 2 », (24) 

where (x 2 ) is related to the average kinetic energy of the b quark inside the B meson, 
see eqs. (56) and (94) in 0, 

(x 2 ) = -^ I (2M B )- 1 (B\bn 2 b\B} = (25) 
3A 3A 



Comparing eqs. ( f23"|) and ([24] ) we conclude that 

p 2 F (2 + p)= A 2 (l + (x 2 )) = ^W 2 (^)(l + (x 2 )). (26) 

Assuming that (x 2 ) is known from a QCD analysis - and to a certain extent this is 
indeed the case - we can use eq. (^6|) to fix the value of p which then may serve as 
an input in eq. (|20"D allowing one to determine pf. Moreover, rewriting eq. fl26|) in 
the form 

l + (x 2 )=7r(2 + p)p- 2 e-^ r 2 (^) (27) 

it is easy to see that the right-hand side is a monotonously decreasing function of p 
having its 111 HX11111 1111 ( ~ 1.6) at p = 0. Thus, in the AC 2 M 2 model an upper bound 
emerges, 

(x 2 ) < | - 1 ~ 0.57 (28) 

On the other hand, the value of /i 2 has been estimated [12] from QCD sum rules 
to be p 2 ~ 0.6 GeV 2 and A ~ 0.4 -h 0.6 GeV. Taking these estimates at face value one 
would conclude that (x 2 ) ~ 0.5 -j- 1. This would mean that the AC 2 M 2 ansatz can be 
made consistent with the real QCD description only for small values of p, i.e. when 
the spectator is relativistic. Further numerical estimates will be presented below. 
The higher moments are obtained in a straightforward way, 

^) ACM = ^p^SK^{^ . (29) 



For even values of n the McDonald function in the right-hand side reduces to an 
elementary one, 

For any value of n the limits of small and large p are 

f 2>£r(^)/v^F atp^O 

\ e ) — * \ n n/2/1 i n(n+2)\ n . n(n+2) P% \ , ■ (31) 

> PFP n/2 (i + -4r J ) = m sv ( 1 + -4 ^) atp- - 



4p I ""spK-^ 1 4 



sp 
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The moments (e n ) given in eq. (|29| ) are related to those of an 'equivalent' AC 2 M 2 
distribution function via the following expression 

(e n ) = A n ((l-x) n > . (32) 

In particular combining eq. (^) for n = 3 and eq. fl2"T| ) we can express the third 
moment, (— x 3 ), via another monotonously decreasing function of p implying that 

< -(x 3 ) < 2 - | ~ 0.43 ; (33) 



for small p the upper bounds of eqs. fl28|), ([33]) actually become approximate equalities. 
One should note that the QCD estimate M (see also 



(-a; 3 ) « 0.03 GeV 3 /A 3 ~ 0.3 0.5 , 
is then quite compatible with the AC 2 M 2 value. 

4. The equivalent 'light cone' distribution function F(x) - the Roman function - 



can be read off from eq. (15) taking into account the definition of the scaling variable 
e = A(l - x), 



n , N 1 A 1 
F R om{x) = —= — exp < -- 
V^Pf 4 



Pf p 




(34) 
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This expression represents how the AC 2 M 2 ansatz can be reproduced by a QCD light 
cone function. The distribution given by the Roman function of eq. fl3~4|) satisfies 
all requirements necessary for the generic distribution function governing inclusive 
decays into massless quarks. As was mentioned, it does not lead to any spectrum 
beyond the kinematical boundary Mb/2 and it exponentially decreases towards large 
negative x. Of course the latter property for F(x) must be considered as a goal in the 
mature QCD description incorporating radiative corrections rather than an obvious 
property; on the other hand it is most natural to require it from the model function 
introduced to describe non-perturbative effects. 

If m sp ^> Pf the Roman function F Rom exhibits a narrow peak around x ~ 
and represents the situation when even the spectator is nonrelativistic. The opposite 
limiting case is more relevant, pp > rn sp . In this case F Rom is rather broad - its width 
is ~ A - and very asymmetric, see Fig. 1. 

5. Next we extend our analysis to inclusive semi-leptonic B decays. The question 
which we want to address here is supplementary to the one treated in Refs. |], 
[TO]] . Namely, we are interested in studying the scaling features of the differential 
distribution 

dT{B -> IviX^/dE^dqo 



8 



where Ei is the charged lepton energy and = (pi+Pu)^ is the total four-momentum 
of the lepton pair. The analysis of Ref. |2| implies that 

dT(B -» W t X u ) _ 2 12( 9o -ff,)(2m^-g 2 ) 

b A TP Aril An ~~ ° ^ 1 > X ^ > 7Z, ~\ ^ 

ahjiaq aqo A [rrib ~ Qo) 

where 

r (6^^ M ) = |K fc | 2 ^|. (36) 

In other words, we predict that the differential distribution - which generically is a 
function of two independent variables, q and q 2 - actually depends only on a single 
scaling combination, 

x = + (37) 
2A(m b - g ) 

apart from a simple kinematical factor. Like in deep inelastic scattering, this scaling 
holds only when both perturbative and higher twist effects are neglected. Both effects 
introduce corrections to the scaling regime which are not discussed here. 

We would like to check whether this scaling feature - the dependence on a single 
scaling variable - persists in the AC 2 M 2 model. To this end we just repeat the steps 
carried out above for the radiative transition, see eq. flT2"|). After averaging over the 
direction of p, the primordial momentum of the b quark, we get for d 3 T/dEidq 2 dqo an 
integral representation quite similar to eq. fll2l) . What is different is the lower limit 
p\ of integration over p 2 which now takes the form 

(38) 




Z=Vb M *V„ 7 ■ (39) 



M 2 +q 2 - 2q M B 
M 2 -q 2 

This expression is sufficient to see that the dependence on q 2 and go enters only via 
the variable e. 

Moreover, e is directly expressed in terms of the QCD variable x of eq. (P7|). 



e = A(l-x) 



l + 0(—(l-x)) 



(40) 



The AC 2 M 2 model, thus, reproduces the correct QCD scaling taking place in the 
leading approximation. 

6. There is, however, an important feature in the true QCD description that is 
beyond the scope of the AC 2 M 2 ansatz: this model, as it is formulated in [||], is 
supposed to be universally applicable to the case of massless and of massive quarks 
in the final state, i.e. to b — ► u and to b — > c transitions, respectively. On the other 
hand, we have shown that the genuine distribution function depends in an essential 
way on the ratio of the quark masses in the initial and the final state, 

F = F(x,rf) , 7 = m 2 /m 2 b . 
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The function is strikingly different in the two extreme cases when 7 <C 1 and 1— 7 <C f ; 
the latter represents the so-called small velocity (SV) limit [|I^j. Let us demonstrate 
that in this limit the AC 2 M 2 model would lead to predictions which differ essentially 
from the true QCD result. 

In the SV limit the velocity of the final quark q is small, \v\ « 1; this can happen 
both in radiative transitions and in semi-leptonic decays; in the latter case it requires 

< (m Q - m q ) 2 - q 2 < m q m Q . 

and in the former 

Am = mq — m q « niq. 

When one retains only terms 0(v °) , the physical spectrum in the B meson decay is 
exactly the same as it was at the free quark level - the 5-function peak resides at the 
same place, and the inclusive probability is completely saturated by one heavy meson 



in the final state ||15|| . Notice that there is no 1 Mb — nib window' - no shift is present 
between the maximal allowed energies of the photon (or lepton) at the quark and the 
hadron levels. 

Modifications of this perfect quark-hadron duality start at the level of 0{y 2 ) |L6|] . 
If 0{y 2 ) effects are considered the height of the elastic peak is changed, and a comb of 
inelastic peaks appears, the height of the latter being proportional to v 2 . This comb 
will lie at E < E Q and will be stretched over an energy interval of order A. ( A is not 
a very relevant parameter in this limit; we can continue to use it, however, just as a 
typical hadronic scale. One could certainly choose another definition of the typical 
hadronic scale.) The integral over the inelastic peaks must compensate the distortion 
of the elastic one - the so-called Bjorken sum rule ||16|| . In the recent paper || we 



demonstrated how this apparently different description emerges from the same QCD 
analysis that leads to the AC 2 M 2 -like prediction for massless final-state quarks. 

Let us compare this picture to the one from the AC 2 M 2 model. Doppler smearing 
yields again one smooth peak whose width is now of order At>. The AC 2 M 2 prescrip- 
tion thus produces nothing that resembles the two-component picture outlined above. 
The failure of this description has a clear origin - it incorporates Fermi motion in the 
initial state but disregards it in the final state where it is now crucial. 

For the same reason in the actual applications of the model in b — > c transitions 
one has to introduce an ad hoc upper cut off for the heavy quark momentum, as was 
mentioned in the Introduction. This cut off is determined by the physical kinematic 
boundary and is, thus, non- universal, i.e. depends on the final quark mass. 

The energy spectrum in the SV limit of QCD is given by a 'temporal' distribution 
function having in general a discrete support representing the higher excited states 
and supplemented by continuum contributions with at least two extra pions in the 
final state. One then has an Isgur-Wise -like |17| description of the spectrum where 



instead of relying on the QCD expressions for inclusive widths one is suggested to 
sum explicitly over the particular final state hadrons. To leading order in \v\ 2 one 
then has only D* (D) and in the next order adds higher excitations Q. Adopting this 
description one then can fit the spectrum by a different phenomenological function not 

2 This is similar to a modification of the original I-W ansatz, incorporating D** states, as now used in 
experimental analyses of semileptonic spectra. 
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directly related to the Fermi motion; it is important however that in the framework 
of the QCD expansion there are certain relations between the moments of the Fermi 
motion distribution function and the parameters of the temporal one (cf. Ref. [[|). 
These relations must be observed in comparing decay spectra with massless and with 
heavy quarks in the final state. 

The mass of the actual c quark is such that one might think at first sight it could 
be safely treated as massless since m 2 /m 2 ~ 0.1. This naive expectation seems to 
be erroneous, though. For the true parameter that enters in semi-leptonic decays is 
rather given by the ratio 

4m 2 m 2 
(m 2 + m 2 — q 2 ) 2 

and its value is typically not very small even in the upper end of the charged lepton 
energy once it is integrated over the invariant mass squared of the lepton pair q 2 . 
The conjecture that the c quark actually lies rather close to the SV limit is supported 
by the observation that the hadronic final state in B meson semi-leptonic decays 
consists, to large degree, of only a D or a D*. (In the exact limit, semileptonic B 
decays would produce only a D or a D* [0.) If so, it is clear that the attempts 
to fit the parameters of the AC 2 M 2 model directly from the data on the inclusive 
b — > c transitions, a quite popular procedure, are not very meaningful from a truly 
theoretical perspective. One can further get convinced that this is indeed the case by 
estimating the value of m c as it emerges from these fits - it turns out to be too high 
from any standpoint. 

It has been claimed in Ref. || that the same description in terms of shape functions 
holds for both b — > u and b —>■ c semileptonic transitions, and that the same hadronic 
matrix elements define it for an arbitrary ratio m 2 /m 2 . It is clear from the analysis 
of Ref. |@ that such claims are erroneous. 

7. The intrinsic inadequacy of the AC 2 M 2 ansatz for describing semileptonic 
b — > c transitions in the SV limit gets obscured once one integrates over the neutrino 
momentum to obtain the charged lepton energy spectrum. Fitting B — > / + X c by 
the AC 2 M 2 model one has commonly set m sp to 150 MeV in a more or less ad-hoc 
fashion; the charm mass together with pp are then used as free fit parameters. In an 
attempt to give this model a somewhat closer connection to QCD one should actually 
adopt a different strategy: namely to require that the difference between the AC 2 M 2 
average of the b quark mass and the charm quark mass satisfy the same relationship 
with heavy flavor hadron masses that the heavy quark expansion yields in QCD for 
m b — m c , i.e. 

(m b ) A CM-m c ~ + (tt) 2m m (41) 

Instead one allows m sp together with pp to float in the fit. It is quite conceivable 
that such an approach would yield values for the fit parameters in better agreement 
with QCD expectations. We do not have at hand the results of such an analysis; still, 
to have an example of the possible distribution function one can use the existing fit 
parameters from CLEO data ||18| , namely pp ~ 282 MeV and m sp ~ 150 MeV. From 
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eqs. (|Sg-gT], H, 0) we get 

A ~ 360 MeV , (vf 2 ) ~ 0.16 GeV 2 

a 2 ~ 0.424 , a 3 ~ -0.28 , a 4 ^ 0.71 , a 5 ~ -1.2 (42) 

and the distribution function itself in terms of the photon energy is shown in Fig. 2. 
Keep in mind that this curve should be taken only as an illustration since the pertur- 
bative gluon corrections, very important for the numerical analysis, are not included 
here. 

One should note that the value of the kinetic energy operator is slightly below the 
lower bound derived in Ref. [Q]. This points again to the intrinsic inconsistency of 
the literal AC 2 M 2 approach with semileptonic decays into charm. The third moment 
of the distribution functions is somewhat less than the estimate of Ref. PJ as well. 
All these facts seem to be correlated with the too high mass of the c quark obtained 
in the AC 2 M 2 fit of such decays. In reality one can obtain better information on the 
distribution function F(x) at small 7 studying the double differential distributions 
and in particular in the region of small q 2 |19|, ||. It is worth reminding however 
that according to the analysis of Ref. Q the calculable 1/mg corrections are really 
important for the fit and they must be taken into account. 

8. So far we have disregarded gluon radiative corrections that actually modify 
the end point shape in an essential way; their interplay with the effects of the Fermi 
motion has been discussed in ref. 0. In radiative decays of the type b — > s+7 the peak 
in the photon energy is still manifest although its height is significantly lowered by 
gluon bremsstrahlung. Both perturbative and non-perturbative effects smearing the 
original narrow line can be described with sufficient practical accuracy by convoluting 
the radiatively corrected parton spectrum with the distribution function describing 
Fermi motion 0]; the radiative corrections must be considered in higher orders of 
perturbative expansion as well, at least in the double logarithmic approximation. A 
more precise treatment including subleading perturbative corrections requires further 
theoretical consideration. We note that it is just this theoretical prescription that has 
been employed in the existing phenomenological analysis (for radiative decays, see 
Ref. [p0||.). In semileptonic decays integrating over the neutrino energy smooths out 
the spectrum and the impact of radiative corrections on the shape is less pronounced. 

The conjecture in Refs. 0, |J that radiative corrections produce relatively small, 
namely 0(a s (rrib)), modifications to the moments of the end-point distribution func- 
tions is certainly incorrect. 

9. In the conclusion of this paper we would like to add some comments on existing 
comparisons of the predictions from the AC 2 M 2 model with QCD . The most obvious 
remark follows from the fact that there are no corrections of order l/m& to total widths, 
as stated first in Refs. [0, §; thus, for example, differences in the lifetimes of B and 
A;, arise first on the 1/m 2 level; the corrections due to Fermi motion are quadratic 
in 1/rrib as well. This does not hold automatically in AC 2 M 2 ; yet by defining an 
effective b quark mass m^ CM as the average over the floating mass of eq. (1) one 
eliminates l/m*, corrections to both the total width and to the regular part of the 
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lepton spectrum. This observation has been noted in a few recent papers |9|, [H 
The above analysis shows that the similarity extends also to the subtle problem of 
the endpoint description for decays into light quarks - contrary to superficial claims 
in Refs. || |J. As a matter of fact the corresponding considerations of Refs. || ^ 
do not refer to the AC 2 M 2 model at all since the crucial effect of the floating mass 
m( has been omitted. As a result the kinematical boundaries were violated and 
the numerical values of the moments obtained there are irrelevant. In contrast, the 



analysis of b — > s + 7 in Ref. |2(| has been done properly: e.g. it preserved the 
kinematic constraints. 

Of course essential differences arise at order l/nig] however one should keep in 
mind that the main point of the AC 2 M 2 model was to take care of only leading 1/mg 
effects that are crucial even in beauty particles. 

In our opinion, this is actually not the main difference, probably not even from 
a practical point of view. For it is important to understand that the value of the 
b quark mass mf CM defined in such a way can actually appear to be different for 
different beauty hadrons Hb and even for the different channels (b — > c vs. b — > u or 
b — > s + 7) in the same hadron. The latter concern follows from the fact that the 
spectrum is shaped by different distribution functions depending on the final state 
quark mass; only the second moments of these distribution functions coincide. This 
difference can actually lead to the emergence of different values for m^ CM beyond the 
nonrelativistic approximation pp <C m sp . 

Once semileptonic or radiative decays of beauty baryons are studied, the second 
concern arises: it is more than conceivable that the literal use of the AC 2 M 2 fit for the 
decays of baryons and mesons would yield values of mf CM that differ already in 1 / mq 
terms; following the standard AC 2 M 2 rules for calculating the width one then would 
arrive at different lifetimes already in the leading non-perturbative approximation. 
For the average mass of the heavy quark plays only a marginal role in the model, 
and thus essentially depends on the assumed shape of the distribution. The whole 
distribution function F(x) can be determined from a precise measurement of the end 
point spectrum, in particular in b — > s + 7 decays. This of course would yield the 
true b quark mass. On the other hand lepton spectra in semi-leptonic decays are less 
sensitive to the exact shape and thus effectively rely more on the model assumptions. 
For in general there is no reason to expect that the particular AC 2 M 2 ansatz for 
the distribution function holds with a good accuracy for real boundstates in the 
presence of potential-like binding forces. For example, a possible and very natural 
generalization of the ansatz eq. ([D would be to add some arbitrary "static" binding 
energy V of the order of Aqcp : 



M B -ml + ^/\p\ 2 + m 2 sp + V (43) 

where V may differ from hadron to hadron. This modification cannot be absorbed into 
a redefinition of m sp and pp. No model can actually provide us with an unambiguous 
prescription for determining the impact of V; QCD on the other hand ensures that 
the effect of V on the total width is cancelled by an analogous interaction in the 
final state. Therefore it seems to be advantageous to impose additional constraint on 
the fitted distribution functions ensuring them to yield the same QCD value of the 
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'average' mass m& independently of a particular ansatz used for their shape. 

We think that the possibility to obtain different mass values for the same heavy 
quark is actually the main physical indication of the inconsistency of the standard 
AC 2 M 2 ansatz with QCD, together with a difference in the underlying Fermi motion 
description for different final state quark masses - which in the case of the SV limit 
is in a sense a rather obvious observation. 
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Figure captions 

Fig. 1. The QCD distribution function F(x) from the AC 2 M 2 model (Roman func- 
tion) for p = 1 (solid line), p = 0.2 (dotted line) and p = 5 (dashed line); p = m 2 sp /p F . 

Fig. 2. The distribution over the light cone momentum 7To + n z corresponding to the 
AC 2 M 2 model with p F = 282 MeV and m sp = 150 MeV . 
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